Large Nc confinement and turbulence 
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We suggest that the transition that occurs at large Nc in the eigenvalue distribution of a Wilson 
loop may have a turbulent origin. We arrived at this conclusion by studying the complex-valued 
inviscid Burgers-Hopf equation that corresponds to the Makeenko-Migdal loop equation, and we 
demonstrate the appearance of a shock in the spectral flow of the Wilson loop eigenvalues. This 
picture supplements that of the Durhuus-Olesen transition with a particular realization of disorder. 
The critical behavior at the formation of the shock allows us to infer exponents that have been 
measured recently in lattice simulations by Narayanan and Neuberger in d = 2 and d — 3. Our 
analysis leads us to speculate that the universal behavior observed in these lattice simulations might 
be a generic feature of confinement, also in d = 4 Yang-Mills theory. 



Many efforts continue to be devoted to the study of 
QCD in the hmit of a large number of colors, after the 
initial suggestion by t'Hooft This, in part, is due to 
the general belief that the large Nc limit captures the 
essence of confinement, one of the most elusive of QCD 
properties. At the same time the theory simplifies con- 
siderably in the large Nc limit: fluctuations die out and 
the measure of the partition function becomes localized 
at one particular configuration, making the large Nc limit 
akin to a classical approximation 2]. 

This paper will discuss the large Nc limit of Yang-Mills 
theory in 2 dimensions, but we have good reasons to be- 
lieve that much of our analysis can be extended to higher 
dimensions. There are many equivalent approaches to 
multicolor Yang-Mills in d = 2. For dcfinitcness, we shall 
refer to the known Q formulation in terms of free random 
variables [1, 0] . This translates d = 2 Yang-Mills theory 
onto large Nc matrix model, where the size of the uni- 
tary matrix is identified with the number of colors. More 
specifically, the basic observable that we shall consider is 
the Wilson loop along a (simple) curve C 



W{A) = ( Fexpi 



c 



(1) 



where the averaging is over the Yang-Mills measure, and 
we have made explicit that W depends in fact only on 
the area A enclosed by C, to within a normalization Q. 
The matrix W is unitary, with eigenvalues of the form 
A — exp{i6) that are distributed, in the limit Nc — s- oo, 
according to an average density p{9,A). 

Recently, Narayanan and Neuberger Q studied the be- 
havior of iy(>l) as a function of A in the large Nc limit. 
They observed that for small loops (which probe short 
distance, perturbative phyics), the spectrum does not 
cover the whole unit circle, but exhibits a gap; in con- 
trast, for very large loops (which probe long distance. 



nonperturbative physics) the spectrum covers uniformly 
the unit circle (gapless phase). Since the crossover region 
is becoming infinitely thin [3] in the limit Nc oo, one is 
tempted to try and explain the transition using classical 
concepts only, which is precisely what we aim at in this 
note. In fact this behavior of the spectrum agrees with 
the order(gapped)-disorder(gapless) transition, proposed 
long ago in the context of large Nc two-dimensional Yang- 
Mills theory by Durhuus and Olesen f6| and based on 
the explicit solutions of corresponding Makeenko-Migdal 
equations 1^. Surprisingly, a similar critical behavior 
has been observed also in d = 3 dimensions and conjec- 
tured to hold in d = 4 large Nc Yang-Mills theory 
In this letter we suggest the general mechanism for such 
a transition, by tracing the complex singularities of the 
eigenvalue flow equation, and demonstrating the appear- 
ance of "spectral shocks" that may signal the transition 
to a turbulent state. 

The spectral density is not available in analytic form, 
but the moments 



Wn{A) =< trVK(A)" >= 



dee"''' pie, A) (2) 



are. An explicit, compact form for these moments is given 
in Ref. [3| in terms of an integral representation 

WniA) = i </-^(l + l/z)«exp(-nA(z + l/2)) 
n J Zm 

= Ul,^_,^{nA)eM-nA/2) (3) 

where the representation of Laguerre polynomials, used 
in the second line, allows connection to results known 
already 25 years ago [1, @]. The Durhuus-Olesen transi- 
tion can be seen by studying the asymptotic behavior of 
these Laguerre polynomials, using a saddle point analysis 
of their integral representation 1^, 11 1 . The result is sur- 
prising: for a loop area below the critical value Ac = 4, 
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the moments oscillate and decay like n~^/^, while for 
A > Ac the moments decay exponentially with n, mod- 
ulo similar power behavior. Both regimes are separated 
by double scaling limit. There exists a vast literature 
on this subject [9|. Here we would only hke to stress, 
that the transition is subtle, and popular arguments that 
two-dimensional confinement is trivial, and perturbative 
(linear potential), can easily lead to paradoxes when, e.g., 
instanton effects are not taken into account 

In order to analyze the nature of the transition, we 
consider the following function [3| 



F{9,A) = i e*''G(e'^^) - 



where G{z) is the resolvent 



G{z,A) = I da- ' 



z — e" 



(4) 
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The function A) is analytic in the complex 0-plane, 
with a discontinuity across the real axis proportional to 
the spectral density, p{6,A) = (l/7r)Imi^(6' - A 
simple calculation yields the explicit expressions [6| 

F{e,A) = ^|^i + 2^„(A)e-™^^ 

= da p{a, A) coi (^-^^ , (6) 

whose imaginary part gives the spectral density in the 
form 



(7) 



in agreement with Eq. ^ above (we have used the fact 
that the moments are real and that p{—0) = p{9)). 

It can be shown that the function F obeys the following 
equation d, d, 



dAF + FdeF = 0. 



(8) 



This is the so-called complex Burgers equation in the in- 
viscid limit. This equation is analogous to the real Burg- 
ers equation of fluid dynamics (with A playing the role 
of time, 9 that of a coordinate, and F of a velocity field). 
The complex Hopf-Burgers equation is omnipresent in 
Free Random Variables calculus This, and similar 
integro-differential equations, also appear frequently as 
one dimensional models for quasi-geostrophic equations, 
describing e.g. the dynamics of the mixture of cold and 
hot air and the fronts between them. Another reason why 
the complex Burgers equation is much studied is the fact 
that the structure of quasi-geostrophic equations resem- 
bles 3-d Euler equations, if one substitutes the velocity 



by the vorticity . Here, we shall take advantage of the 
abundant mathematical studies of the complex Burgers 
equation to suggest a connection between signals of tur- 
bulence (commonly associated with the blow-up of the 
solution in finite time) and the problem of spectral flow 
of the eigenvalues of Wilson loop operators. We shall in 
particular adapt the proof of the blow-up from Ref . 1J| , 



using the method of complex characteristics, and trac- 
ing singularities in the complex plane. The observation 
of shock formation in the complex Burgers and similar 
equations was confirmed by other methods, see e.g. [l^. 

In the present case, the method of characteristics pro- 
vides the following implicit solution in terms of an aux- 
iliary function ^(A, 0): 

F{A,9)=Fo{aA0)), 9 = ^ + AFo{0, (9) 
where Fq is given by the initial condition 



Fo{9)^F{A = 0,t 



(10) 



This implicit solution is well defined as long as the map- 
ping between 9 and ^ remains regular. However singular 
points occur when d9/d^ = 0, that is, for solution of 



l + AF^{ic)=0. 



(11) 



A lot of information on the solution, and in particular 
the occurrence of a blow up, can be inferred from the 
study of the location of these singular points. 

We begin our analysis with the ordered (gapped) state, 
and recover known results (see e.g. the second Ref. [lo| ) 
in an easy and direct way. We start from an initial con- 
dition peaked at eigenvalue A = 1, i.e. p{9,0) — S{9) 
(corresponding to A = in Eq. Q). Then, from ©, 
Fo{9) = l/2cot6'/2, and from ^ sin^ = A/4. In 
the vicinity of the singularity the characteristics behave 



^cf^T/A~l 



(12) 



where 0^ = + ~ ^/4) and = 2 arcsin(VA/2). 

The spectral density is easily deduced from the imaginary 
part of F or equivalently Fq (see (ITO|) '). A simple analysis 
then reveals that p{9,A) = when 9 > 9c, while for 
9 ^ 9c, p{9,A) ^ \J9c~ 9. In other words, 9c determines 
the edge of the spectrum which, as long as A < 4 presents 
a gap. 

At the closure of the gap, i.e., for A = A, the second 
derivative in the expansion (jl2p vanishes, and we have 
instead 



(13) 



with the spectral density p{9, A) ~ (tt — 9)^ for 9 ^ 
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From this point on, as A keeps increasing, Sft^c remains 
equal to tt, but the initial singularity splits into two com- 
plex conjugate ones that move away to zbicxi, leading 
eventually to a uniform spectral density at large A. This 
behavior is reminiscent of the turbulent inverse spectral 
cascade (alike in two-dimensional turbulence which, 
as A grows, suppresses higher Fourier modes, leaving in 
the A = oo limit only the longest wavelength mode cor- 
responding here to the constant density p{6) — 1/2tt. 

To demonstrate the inverse spectral cascade, we con- 
sider, following a small perturbation of the form 



p{0, Ao 



2ecos0) +0{e^), 



(14) 



where Aq >> 1. This is of the form ([7]) with wi ~ e, 
and all other moments vanishing exponentially. It fol- 
lows from © that Fo{0 = |(1 + 2eexp(-i^)), 9 ^ C + 
{A — Ao)Fo{^), and a singularity occurs when exp(z^c) = 
—€{A — Aq). In the vicinity of the singularity 



+ ^(e-ec)^ 



(15) 



where 6c ^ £,c + {A — Aq)Fq{^c)- We have two solutions 
depending on whether A > Aq or A < A^. In the first 
case 

n-tll~^^-^ + \ne{A-Ao)], (16) 



whereas in the second case 

• (1 , ^0-^ 
M 1 H ^ 



+ lne(^o - A) 



(17) 



In the first case, the singularity is initially (when A = Aq) 
at 9c = n + ioo, and for small e (e < 1/2), it remains 
complex, and returns to +ioo as ^ — > oo: asymptotically 
the efltect of the perturbation vanishes, and only the con- 
stant mode survives. However in the second case, the 
initial singularity at -|-ioo moves towards the real axis as 
A decreases, and reaches it in a finite time A* given by 



= 1 



An -A* 



+ lne(Ao-^*). 



(18) 



At this point we have a blow-up [1J| , of the type already 
encountered earlier at the closure of the gap. In fact, 
this phenomenon is quite generically associated to the 
motion (as A decreases) of complex conjugate singular- 
ities (square root branch points) towards the real axis, 
that eventually eventually merge on the real axis (into a 
third order branch point). 

Th is g eneric phenomenon has an exact analog in op- 
tics [16| , where light rays play the role of the character- 
istics. The singularities lines are there the caustics, and 
the merging of singularities discussed above corresponds 
to the merging of two "fold" caustics into a "cusp" (in 



the terminology of catastrophe theory). The essence of 
the cusp singularity is captured by the Pearcey's function 

P{Lri)^ dtexp{t{f/4 + ^ty2 + 7jt)) (19) 



whereas the square root type singularity is described by 
the well known Airy function. This analogy can be used 
in order to explain the origin of the universal function 
for Yang-Lee zeroes for d > 2 dimensional Yang-Mills 
obtained in Ref. The explicit comparison will be 
presented elsewhere (23 |. 

The location of the singularities in the complex plane 
also determine the asymptotic behavior of the moments, 
according to general arguments When 9c is close 

to the real axis, one may write 6c = 9*+ iA{A), with 9* 
real; then the singularity is of the form {9 — 6'*)^ , where 
fi is not an integer (branch point singularity). It follows 
that the large n Fourier coefficients, i.e., the moments 
Wn, behave as 



Wn = kr(^+i)e-' 



(20) 



Thus, the position of the singularity determines the width 
of an analytical strip, controlled by the value of A(A). 
In the gapped phase, the singularity is of the square root 
type, but always on the real axis. Then, since A = 0, 
the moments only oscillate, accompanied by power law 
^-3/2 _ -j-j^g gapless phase, the same power law is ac- 
companied by exponential damping, since in this case 
the singularity is complex (with a non vanishing A{A)). 
At the critical point, a generic cubic singularity (see 
Eq. (fT3|l and [ll,[3]) appears in a narrow interval where 
1 << n << 1/A{A) (just approaching the singularity 
from above) and the scaling goes like 



1-4/3 



(21) 



This generic picture agrees with detailed calculations 
using the analytic form of the moments given above, 
based on asymptotics of Laguerre polynomials 0, 113, U | ■ 
We would like to stress however that the method of trac- 
ing complex singularities is general and therefore we ex- 
pect it to hold also in higher dimensions. The "speed" 
with which the singularities move may depend on the di- 
mensionality, and possibly on other features of the loops, 
but for large loops it is natural to expect that the domi- 
nant control parameter will remain the area of the loop. 
We may therefore speculate that, also in higher dimen- 
sions, the disordered phase may be caused by an inverse 
spectral cascade. 

One may also argue that the universal nature of the 
critical behavior, conjectured and observed in Ref. 0], 
may find its origin in the fact, that since the critical be- 
havior happens in a very narrow analyticity strip, one 
can expect universal (ergodic) behavior, alike in several 



4 



models of disorder in mesoscopic physics. If this is the 
case, one may expect that simple schematic models (e.g. 
matrix models) may define this class of universality. This 
seems to be indeed the case. A particular illustration is 
provided by the matrix model proposed by Janik and 
Wieczorek [l^, hereafter JW model. The model stems 
from the general construction of multiplicative free evolu- 
tion [20I , where increments are mutually free in the sense 
of Voiculescu. The unitary realization in JW model cor- 
responds to matrix value unitary random walk, where 
the evolution operator is the ordered string of consecu- 
tive multiplications of infinitely large unitary matrices 



K 



W 



(22) 



where Ui = exp ^^Jt/ KHi , with Hi a hermitian random 
matrix. The model is a random matrix generalization 
of the multiplicative random walk performed in K steps 
during "time" t. In the continuum limit K 00, the 
model is exactly solvable. The solution coincides exactly 
with the two-dimensional QCD, provided one identifies t 
with the area of the Wilson loop, modulo a normalization 
a. 

Recently, Neuberger and Narayanan have observed, 
that large Nc Yang-Mills lattice simulations in c? = 2 and 
d~?> demonstrate the same critical scaling at the closure 
of the gap as in the JW model and have conjectured that 
this model establishes a universality class for d — 4 large 
Nc Yang-Mills as well. That the whole dynamics of com- 
plicated non-perturbative QCD can be reduced in some 
spectral regime to matrix model is not new - a notable 
case is the universal scaling of spectral density of Eu- 
clidean Dirac operator for sufficiently small eigenvalues, 
where the spectrum belongs to broad universality class 
of corresponding chiral models 22]. In fact the present 
analysis leads us to expect that in the very narrow spec- 
tral window around A = — 1 a universal oscillatory regime 
precedes the formation of the spectral shock, in qualita- 
tive analogy to similar spectral oscillations of the quark 
condensate before the spontaneous breakdown of chiral 
symmetry, based on the Banks-Casher relation [2^. 

In this note, we have proposed to view confinement- 
deconfinement transition in multicolor Yang-Mills as a 
order-disorder phenomenon, where the transition to dis- 
order is caused by (inverse) turbulence in the spectral 
flow of Wilson loop operators. This picture corroborates 
the picture of order-disorder transition that Durhuus and 
Olesen have envisioned many years ago, and supplements 
it with a detailed model for building the disorder, based 
on the development of the inverse spectral cascade. It 
will be interesting to incorporate into this picture the 
effects of finite Nc and matter fields. We expect, that 
finite Nc effects will contribute to the appearance of an 
effective spectral viscosity which will smoothen the 
shocks, but will not destroy them 
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